Meissner effect in a bosonic ladder 
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We investigate the effect of a magnetic field on a bosonic ladder. We show that such a system 
leads to the one dimensional equivalent of a vortex lattice in a superconductor. We investigate 
the physical properties of the vortex phase, such as vortex density and vortex correlation functions 
and show that magnetization has plateaus for some commensurate values of the magnetic field. 
The lowest plateau corresponds to a true Meissner to vortex transition at a critical field Hd that 
exists although the system has no long range superconducting order. Implications for experimental 
realizations such as Josephson junction arrays are discussed. 
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The effect of a magnetic field on interacting particles is 
a long standing problem. A spectacular case is provided 
by type 11 superconductor, in which the magnetic field is 
totally expelled below Hc^ , whereas a vortex state exists 
for H > Hc^ . This behavior, however, is obtained from 
the Landau-Ginzburg equation, and it is important to 
know what happens when interactions and fluctuations 
have more drastic effects, such as in one dimensional sys- 
tems. Indeed, in a one dimensional conductor at T = 0, 
although there is no long range order, superconductivity 
in the sense of infinite d.c. conductivity can nevertheless 
be presentEJ. For a one dimensional chain, there is no or- 
bital effect, and this question is not relevant. However, 
a system made of a finite number of coupled chains (a 
ladder), is still one dimensional (no long range order can 
exist) but orbital effect of the magnetic field is present, 
opening the possibility of such a transition. 

Beyond its own theoretical interest the investigation 
of the effect of a magnetic field on ladder systems is also 
of direct experimental releYance, due to the various re- 
alizations of such laddergja'cl. Fermionic ladders can be 
superconducting both from attractive (s-wave) and re- 
pulsive interactions (d-wave). In the attractive case, the 
system is close to standard superconductors where pairs 
of fermions can hop from one chain to the other leading 
to a Josephson coupling, provided the applied magnetic 
field is smaller that the spin gap. The system can thus 
be dewajibed as a bosonic ladder. Josephson junction 
array^pQ provide-.aJso a very direct realization of such 
a bosonic laddeiuEl and are thus the prime candidates 
to observe these effects. This problem of Josephsou. lad- 
ders has been investigated previo us I j^f in the classicaB and. 
quantum limit both analyticallyll30 and numericallyEHI 
in the high field limit of half a flux quantum per pla- 
quette for one dimensional situatixmSL. For the Josephson 
two leg ladder, it has been shownEElEJ that a true transi- 
tion exists between a commensurate and a vortex phase, 
however the detailed behavior of the vortex phase and the 
effects of commensurability of the magnetic field remain 
to be understood. 



In this paper we investigate the effect of a magnetic 
field directly on the bosonic two leg ladder. We study 
the Meissner-vortex transition in this system, and the 
nature of the vortex phase. We show that plateaus in the 
magnetization for commensurate values of the field exist, 
in a way similar to the Mott transition in one dimension 
for commensurate values of the filing. We analyse the 
consequences for transport in such bosonic ladders. 

The lattice Hamiltonian of the bosonic two leg ladder 
in a magnetic field is: 

i,p=l,2 

i 

+ U^Ui^pirii^p - 1) + Vni^in^^2 (1) 

where the density rii^p = b] ^bi^p and the magnetic field 
is introduced via the Peierls substitution. In addition to 
their parallel ty and perpendicular t± hopping the bosons 
repel via an on-site U and an inter-chain V repulsion. 
The presence of a uniform magnetic field implies that: 

A - dl^ Aj_{i + 1) + a^||.i(i) - A±{i) - a^||.2(i) = ^ 

(2) 

Where $ is the flux of the magnetic held through a 
plaquette. r^he low energy properties of (|l]) are more 
transparentll3 using a "bosonized" representation of the 
boson operators^. One obtains: 

^ = E [ ^[uKinHp^ e*A^^,pr + ^{d^^p)'' 

— dx cos{9i - 0-^ + e* Aj_{x)) + 
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where (pi and Trlli = V^^i are conjugate variables and 
the boson annihilation operator is given by: 



il}p{x ~ no) = — ^ 
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(4) 



where a is the lattice spacing along the legs of the ladder. 
u is the sound velocity of the collective density oscilla- 
tions. The linear nature of the spectrum is the signature 
that true superfluiditw exist in a single chain even with- 
out long range ordem. K is the Luttinger parameter, 
directly related to the compressibility of the system, and 
incorporating all microscopic interaction effects. For a 
simple onsite repulsion 1 < K < oo, K = oo being free 
bosons and K = 1 hard core bosons. For more gen- 
eral (longer range) interactions all values of K are in 
principle allowed. If one introduces the symmetric and 
antisymmetric combinations ^^.q = '^^^^ (and similar 

combinations for 11 and A), Eq. (||) readsEl: 
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As describes an Ahronov-Bohm Flux threading the sys- 
tem. In the following, we assume that there is no such 
flux, so As = 0. We also assume that Ka > 1, so 
that the term cos\/80a is irrelevant. This leads to an 
Hamiltanian similar to the one derived for the Josephson 
laddeitjlid. The total and antisymmetric longitudinal 
current ^ = ji ± j2: and the current perpendicular to 
the ladder are given by (i/ = s, a): 



A, 
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Hg gives a simple dynamics, the behavior of Ha is 
richer. In the absence of a magnetic field. Ha is a sine- 
Gordon Hamiltonian. For i^T > 1/4, it develops an order 
in the field da with (da) ~ and a gap in the excitation 
spectrum (^L^y2-i/i2K) ^ rpj^j^ describes the lock- 

ing of the relative phases by Josephson coupling. Let us 
now consider the effect of a magnetic field such that the 
flux per plaquette $ = % 2 ^ where p, q are two mutually 
prime integers. Using the gauge Aa = 0, A± ~ we 
can write: 
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FIG. 1: The phases of the bosonic two leg ladder in the pres- 
ence of a magnetic field, (a) The low field phases. Current 
only exist along the legs of the ladder leading to screening of 
the appfied magnetic field, (b) the high field phase. Current 
exist both on the legs and the rungs of the ladder leading to a 
vortex lattice (here shown with a vortex every three sites). At 
commensurable fiux the vortex lattice is perfectly ordered (at 
T = 0) whereas away from commensurability the vortex sys- 
tem has only quasi long range positional order and no perfect 
crystalline order. 



The Hamiltonian (^ is similar to-the one describing the 
Mott transition in one dimensionEj, where any commen- 
surability can in principle give a transition depending on 
K. Indeed, even if the presence of an oscillating phase in 
(^ makes tj_ naiveW irrelevant, perturbation theory to 
order q in showsE3E3 that a term: 



H„ = ^ I dx cos qV2da 
na 



(10) 



exists with gqa/{'KUa) ^ {t^a/nUaY in perturbation. 
More generally this term is allowed by symmetry and 
is present in the Hamiltonian even beyond perturba- 
tion theory. The term (|lO|) opens a gap ^p/q — 



t± c 



" ' provided Ka > \. The expectation 

value of Oa in the gapped phase is {9a) — {2k + l)'K / {qy/2) 
where k = 0, . . . , g — 1. This corresponds to a g— fold 
degenerate ground state resulting from the breaking of 
the discrete tseHislation symmet ry u ii p direct analogy with 
Mott systemsEJ or spin systemsllSEj. The corresponding 
expectation value of j_L (x) is: 
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giving a periodic pattern of the transverse currents of 
period qa as shown in Fig. |^. This pattern corresponds 
to a vortex lattice phase pinned to the microscopic lattice 
with p vortices in a supercell of q sites, leading to an 
average vortex density = ^ . Let us now take a fiux per 
plaquette which is an irrational multiple of the quantum 
of flux. If this value is close to a rational number with 
a small enough denominatorwe can decompose e*^ = 
27r- -I- 5$. The most convenient gauge choice is then 

A± = {2TTe*p)/q and Aa = (5$)/(a\/2). Introducing the 
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field Pa conjugate to 9a, irPa = dx(t>. Ha becomes: 
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The deviation of the flux per plaquette from a ratio- 
nal value in the unit of the quantum of . fl i -ix. ^causes 
a commensurate-incommensurate transitioriB3EJ'c3: At 
low deviation, the system is in its commensurate phase, 

^{Oa) = (2fc + l)7r/g. This implies (j,) = -uK^-^^ 
and no modification of the perpendicular current {5i±) = 
and {5jA^{x)5jj^{x')) ~ e-l=^-^'l/«° with = u/A^. In 
particular, for low magnetic field (p = 0, g = 1), there 
is a Meissner current circulating only on the edges of 
the system (seepE'ig. |]) and proportional to the apphed 
vector potentially. For a large enough deviation from ra- 
tionality, ^(5$ > Ap/q, there is a transition to a phase 
where the cosgv^^a operator is irrelevant. In this phase 
{dxOa) / {t^V^) is non zero an defines an order parameter. 
In the small applied field casellJ, this order parameter is 
just the vortex density pv- In the general case, this order 
parameter measures the deviation Spy = Py ~ | of the 
vortex density from a rational value. The average trans- 
verse current is zero in this phase. However, transverse 
currents correlations have an oscillating component that 
decays with a power law as 



cos[27rpy (a 
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where K* is the renormalized Luttinger parameter. It 
would be, worthwhile to check in a numerical simulation 
such asli3 whether such current correlations are present. 
There is a pattern of transverse currents alternating along 
the rungs of the ladder. One can identify the solitons in 
9a with vortices surrounded by circulating currents (see 
figure |l]). If a vortex is at position x, it causes a jump of 
2tt in 9 1 — 6*2 at this point. This enables us-to write the 
density of these one-dimensional objects aalj: 
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When the density of vortices is p/q, the potential energy 
term ( |To| ) can be interpreted as the coupling of the vortex 
density to a pinning potential of period q lattice spacings. 
Therefore we recover the interpretation in term of pin- 
ning of the vortices by the underlying microscopic lattice 
and the analogy with the Mott transitionta in d = 1. 
Eq. ([l^ ) gives as vortex current jv{x) — The ap- 

plication of a voltage Vi — V2 between the chains gives: 
{dx(ba) = e*^^^ and - F2 = E^b = jy, which 
is just Faraday's law. This completes the identification of 
the incommensurate phase as an unpinned vortex lattice. 
Because of the linearization of the spectrum in (||) there 
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FIG. 2; Sketch of the staircase in the magnetization of the 
bosonic or Josephson ladder. On the figure, Ka = 4, so that 
only the plateaus obtained for pjq — 0, 1/3, 1/2, 2/3, 1 survive 
to the quantum fluctuations. Note that the width of the 1/3 
and 2/3 plateaus is already extremely reduced compared with 
the width of the 1/2 plateau. 



is no Laplace force on vortices. This pathology could be 
cured by putting bacLthe band curvature. 

In the classical vortex lattice phase is obtained 

each time the flux per plaquette in a rational multiple 
of the quantum of flux, leading to a devil's staircase 
structure in the behavior of the magnetization. Here, 
the quantum fluctuations wipe out the large fractions 
for which q^ > AKa, so only some plateaux remain as 
shown on Fig. ^ As quantum fluctuations are reduced, 
more and more plateaus in the magnetization curve are 
formed. It is interesting to study how the width of a 
magnetization plateau increases as a function of Ka or 
tj^. Standard RG calculation of the gap show that the 
width of a plateau behaves as exp(— C/ yjt± 
to the threshold. For <I> 



i^) close 

^ , a vortex lattice phase has 
been obtained in numerical simulationsEHI. It remains to 
be seen whether such phase is a Luttinger Liquid or a 
pinned vortex lattice. 

Close to the transition the coefficient K* 



IS univer- 



sal and Ka = 5^/2. It is-|pQssible in the limit 5$ 



5^ci{p/q) to rcfcrmionizi 
tex density {pv) = 



yielding an average vor- 



and vortex current 
The transition from the 



pinned vortex lattice state to the vortex liquid state is 
thus continuous. The behavior of the average current as 
a function of flux is plotted on figure ^. This is reminis- 
cent of the behavior of a Type II superconductor, howevei; 
the growth of the magnetization is here quite differentE2l 
from the standard B = H + M ^ \n^{H-H 1) 

How could the Meissner and the vortex lattice phase 
be detected in an experiment ? First, hc^ = -fg- must be 
smaller than the spin gap. Otherwise, the magnetic field 
would break the Cooper pairs before vortices are nucle- 
ated, a situation reminiscent of type I superconductivity. 
Second, in a standard material, a ~ lAso that /ici ~ 
75-7-5- ~ 2.1 X lO^T making the vortex lattice unobserv- 
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FIG. 3: The behavior of the screening current with "!>. For 
$ < $c, the system is in the Meissner phase and the screening 
current increases linearly with "1>. For "1> > $c, the screening 
current decreases with the applied field. 

able. For Josephson junction ladders, a ~ 1/xm, giving 
a typical /id — 2.1 x 10^'^T which is easily achieved. In 
principle, the formation of the vortex lattice in this sys- 
tem could be observed by magnetization measurements 
or by resolving indiiiLdual vortices as was done for other 
mesoscopic systemsEj. The Meissner transition also af- 
fects transport properties at a commensurate filling of an 
integer number of Cooper pairs per island. In this case, 
a cos 20s cos 20a term is present and favors an insulat- 



ing stateEil. The interchain Josephson coupling cos \/2^a 
competes with this tepaa making the system conduct- 
ing when large enoughca. The application oi H > Hd 
suppresses the interchain Josephson coupling turning the 
ladder insulating in the vortex phase. The Meissner tran- 
sition is thus accompanied by a superconductor-insulator 
transition in this case. At incommensurate filling, the 
magnetic field still affects transport properties if there 
is an artificial defect in the ladder, which adds a term 
VicosV2(j)s{x = Q)cos^/2(j)a{x = 0)-hV'2Cos V80s(a; = 0) 
to the Hamiltonian. In the Meissner phase, the field 
cos \/2(j)a is disordered and the term dominates trans- 
port, li K > 1 this term is irrelevant and the conduc- 
tance tends at low temperatures to the quantum of con- 
ductance. In the vortex phase on the other hand, the 
field cos \/20a has quasi long range order and the Vi term 
dominates transport, leading (for X < 2) to an insulat- 
ing state at zero temperature. More generally, even when 
K < \ OY K > 2 the conductance has a qualitatively dif- 
ferent temperature dependence in the Meissner and in 
the vortex phase. 
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